Small-Network Approximations for Geometrically Frustrated Ising Systems 
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The study of frustrated spin systems often requires time-consuming numerical simulations. As 
the simplest approach, the classical Ising model is often used to investigate the thermodynamic 
behavior of such systems. Exploiting the small correlation lengths in frustrated Ising systems, we 
develop a method for obtaining a first approximation to the energetic properties of frustrated two- 
dimensional Ising systems using small networks of less than 30 spins. These small networks allow 
much faster numerical simulations, and more importantly, analytical calculation of the properties 
from the partition function is possible. We choose Ising systems on the triangular lattice, the 
Kagome lattice, and the triangular Kagome lattice as prototype systems and find small systems 
that can serve as good approximations to these prototype systems. We also develop criteria for 
constructing small networks to approximate general two-dimensional frustrated Ising systems. This 
method of using small networks provides a novel and efficient way to obtain a first approximation 
to the properties of frustrated spin systems. 

PACS numbers: 75.10.Hk, 75.10.-b, 75.10.Jm 
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I. INTRODUCTION 

The phenomenon of frustration in condensed-matter 
systems has been studied for over 50 years. In 1936, 
Pauling first noted that frustration is present in the struc- 
ture of ice. The tetrahedral structure of ice allows mul- 
tiple possible locations for the hydrogen atoms, giving 
rise to about (3/2)^ ground-state configurations for a 
total of N water molecules.^ Pauling's prediction was 
later confirme d by the experimental work by Giauque 
and co-workers .l^^^ Since then, frustration in condensed- 
matter systems has been an active area of research. Early 
theoretical investigations include the study of the anti- 
ferromagnetic triangular Ising lattice by Wannier, who 
showed that the system is disordered at all temperatures 
and the ground state entropy is 0.323A^/c, where N is the 
number of spins and k is the Boltzmann constant.^ In 
recent years, new frustrated materials have been discov- 
ered and characterized following advancements in fabri- 
cation and measurement techniques. Most notably, the 
ground state spin configurations of rare earth pyrochlores 
Ho2Ti207 and Dy2Ti207 have been found to have a one- 
to-one correspondence to the structural configurations of 
ice, and thus they are aptly described as the "spin ice" W^ 
Adding to their interest, Castelnovo et. al. recently pro- 
posed that magnetic monopoles emerge in the frustrated 
spin ice system,^ and this has been supported by a few re- 
cent experiments.^^H^ Apart from three-dimensional frus- 
trated systems such as spin ice, two-dimensional frus- 
trated systems with a range of different geometries can 
now be artificially constructed with nanometer-size mag- 
netic islands^^ or closely-packed colloidal spheres. ^^ In 
addition, it has been proposed that frustrated spin sys- 
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tems can be constructed by trapping cold atomJHI or po- 
lar molecules^^ in an optical lattice. An important leap 
towards the experimental realization of such systems has 
been achieved by Simon et. al. by succesffully demonstrat- 
ing antiferromagnetic spin chains in an optical lattice.^ 

Frustrated magnetic systems exhibit many interesting 
properties. Most notably, they have multiple degener- 
ate ground states, which give rise to a non-zero entropy 
at absolute zero temperature, violating the third law of 
thermodynamics.^^ In addition, even in the regime when 
kT is much less than the energy scale of the spin-spin 
interactions, there can still be significant fluctuations 
within the system. ^^ Moreover, many frustrated systems 
possess rich phase diagrams as temperature or exter- 
nal magnetic field strength is varied. The phases dis- 
played in these sytems can be magnetically ordered, par- 
tially ordered, or completely disordered.^ In addition to 
their theoretical interest, frustrated magnetic materials 
may have novel technological applications from micro- 
electronics to drug delivery. For example, efficient and 
environmentally friendly magnetic refrigerators may be 
constructed with frustrated systems using the technique 
of adiabatic demagnetization^^ and new technologies in 
advanced magnetic-recording devices may also be built 
with frustrated materials .^^ A better understanding of 
the behaviors of frustrated systems may even allow us to 
gain insights in fields beyond condensed matter physics. 
For instance, it has been suggested that the folding of a 
protein into a biologically funtionable structure is a re- 
sult of the natural ability for the protein to resolve the 
frustrated couplings. ^^ 

As the simplest approach, the thermodynamic proper- 
ties of geometrically frustrated systems can be studied 
with the Ising model. For systems in which classical spin 
fluctuations dominate, the Ising system can be a very ac- 
curate physical model. For other systems in which quan- 
tum fluctuations dominate, the Ising system becomes less 
accurate, ^^ but we may still use it to obtain a first ap- 



proximation. Currently, there are two primary methods 
for studying frustrated Ising systems: exact analytical 
methods applied to infinite systems and Monte Carlo sim- 
ulations applied on very large systems and extrapolated 
to the infinite-system limit. However, exact analytical so- 
lutions are not always possible, while Monte Carlo simu- 
lations can often be very time-consuming. Therefore, we 
are motivated to develop an efficient first approximation 
method for geometrically frustrated Ising systems. 

One of the key features of geometrically frustrated sys- 
tems is their small spin-spin correlation lengths. In the 
absence of long-range correlations, many properties of the 
system are determined by the local geometric network in 
the vicinity of each spin. By arranging a small number of 
Ising spins in an appropriate network, we find that the en- 
ergetic properties of the extended two-dimensional Ising 
systems can be reproduced with surprising accuracy. Not 
surprisingly, the shorter the correlation length in the ex- 
tended system, the better it can be approximated by a 
small network. There are two main advantages of us- 
ing such small- network approximations. Firstly, since 
the small networks that we develop have less than 30 
spins in general, it is much more efficient to do Monte 
Carlo simulations on the small networks. Secondly, for 
networks of less than about 25 spins, it is possible to cal- 
culate the thermodynamic properties of the system di- 
rectly by utilizing the Boltzmann distribution with mini- 
mal computational power. This method of small-network 
approximation may provide us with a new and efficient 
first approximation to the properties of frustrated sys- 
tems. When an appropriate small network is used, both 
exact analytical methods and Monte Carlo simulations 
are much simpler than for the associated extended Ising 
systems. 

In this paper, we present an approximation technique 
for frustrated Ising systems; our goal is to find small 
Ising networks with less than 30 spins that accurately 
approximate the energetic properties of extended Ising 
systems in the thermodynamic (infinite-size) limit. We 
illustrate this technique on the two-dimensional triangu- 
lar, Kagome, and triangular Kagome lattices. Since the 
energetic properties of these systems can be well repre- 
sented by the specific heat vs. kT profiles, we attempt 
to find small networks that accurately approximate the 
specific heat of the extended systems and develop general 
criteria for constructing good small-network approxima- 
tions. The paper is organized as follows. In Section [TT| 
and |III[ we describe our models and the numerical tech- 
niques that we use to study the models. In Section |IV[ 
we compare the specific heat of particular small networks 
to that of the extended Ising systems on the triangular 
Kagome lattice, the Kagome lattice and the triangular 
lattice. In Section [V| we elucidate the general criteria 
for constructing a small network that well- approximates 
the specific heat of the extended two-dimensional lattices, 
and in Section |Vl| we discuss the advantages for using 
small-network approximations over carrying out simula- 
tions on a small piece of the extended lattice with pe- 



riodic boundary conditions. Finally, in Section |VTTj we 
offer some concluding remarks. 



II. MODEL DESCRIPTION 

A. The Extended Ising Models 

In this paper, we discuss the possibility of using small 
Ising systems to approximate the specific heat of geo- 
metrically frustrated two-dimensional Ising systems. We 
have chosen the Ising systems on the triangular lattice, 
the Kagome lattice and the triangular Kagome lattice 
as prototypes of extended two-dimensional systems. The 
network structures of these three lattices are schemati- 
cally displayed in Fig. [l] in which a solid or open circle 
represents a site for an Ising spin, and a single line repre- 
sents the nearest-neighbor bond between the spins. We 
note that the triangular Kagome lattice is different from 
the other two lattices in that while all spins on the tri- 
angular lattice and the Kagome lattice are geometrically 
equivalent to one another, there are two kinds of sites 
with distinct geometries on the triangular Kagome lat- 
tice. We have distinguished the two kinds of spins using 
open and solid circles in Fig. [ijc), and we shall refer to 
the spins represented by the open and solid circles as 
a-spins and 6-spins respectively. 




FIG. 1. Schematic representations of the structures of ordi- 
nary two-dimensional Ising systems on (a) the triangular lat- 
tice, (b) the Kagome lattice and (c) the triangular Kagome 
lattice. The solid and open circles represent spin positions 
and the shaded areas represent a unit cell. 



The Hamiltonian for an Ising system is given h^p^ 



The Hamiltonian of a small network is given by 



{ij) 



(1) 



where the indices i and j label the lattice sites, ai = ±1 
is the spin variable at site z, representing the up and down 
directions of the Ising spin, and Jij is the interaction be- 
tween the spins on i and j sites. A positive value for Jij 
represents antiferromagnetic interactions while negative 
Jij represents ferromagnetic interactions. The summa- 
tion runs over all nearest-neighbor pairs (ij). 

In the case of the triangular lattice and the Kagome 
lattice, all spins are geometrically equivalent and all the 
nearest-neighbor interactions are of the same magnitude. 
Thus, Jij = J and we can rewrite the Hamiltonian as 



H = jy^dicrj 

{ij) 



(2) 



In the case of the triangular Kagome lattice, there 
are two kinds of nearest-neighbor interactions: the 
aa-interaction and the a6-interaction. We denote the 
strength of the two kinds of interactions as Jaa and Jab 
respectively, and the Hamiltonian of the system is given 

by 



f^ — 'Jaa / ^ ^a,i 
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(3) 



where dx^i denotes the spin with index i on the x-lattice 
{x = a, 6). The first summation runs over all nearest- 
neighbors (ij) within the a-sublattice and the second 
summation runs over all nearest-neighbors {ij) between 
the a- and the 6-sublattices. 



B. The Small Ising Networks 

The small Ising systems that we have constructed to 
approximate the extended lattices in the previous section 
span a range of structures. Each consists of a group of 
less than 30 spins connected in a specific network. In 
many cases, the network of a particular small Ising sys- 
tem can be represented by a polyhedron, on which the 
vertices represent the spin positions and the edges rep- 
resent the nearest-neighbor bonds. We will display the 
structures of the small systems as we discuss the small- 



network approximations in Sec. IV 



As in the extended Ising systems, each spin position on 
the small networks has a spin variable of cr^ = ±1, repre- 
senting the up and down orientations of the spin. Since 
most of our small networks can be represented by three- 
dimensional structures such as polyhedra, the terms up 
and down may cause potential confusion due to the sur- 
face curvature of the structures. In our case, up and 
down represent the absolute up and down directions in 
spin space and not the directions perpendicular to the 
surface of the polyhedron at any particular spin position. 



n = ^JijCJiCTj 

{ij) 



(4) 



where Jij is the interaction between spins on sites i and 
j, and the summation runs over all the bonds in the net- 
work. 



III. NUMERICAL TECHNIQUES 

In this section, we introduce the numerical methods 
involved in our study. Section HI A gives the numeri- 
cal details of the Monte Carlo simulation, which is used 
to calculate the properties of both the extended lattices 



and the small networks. In Section IIIB[ we present the 
method for calculating the specific heat exactly from the 
partition function; this method is computationally inex- 
pensive for systems with less than about 25 spins. In Sec- 
tion [TTrCj we introduce a "deviation index" to quantify 
how well each small network approximates the thermo- 
dynamics of the extended lattice system. 



A. Monte Carlo Simulation 

Frustrated Ising systems exhibit little long-range order, 
and therefore a single-spin-flip algorithm is suitable to 
statistically sample the microstates of these systems. In 
view of this, we employ the Metropolis algorithm to study 
the Ising systems on both the extended two-dimensional 
lattices and the small networks .^^^^ For ordinary ex- 
tended two-dimensional Ising systems, we carry out the 
simulations using systems of L x L unit cells under pe- 
riodic boundary conditions. We use L = 30, 18 and 10 
for the triangular lattice, the Kagome lattice, and the 
triangular Kagome lattice, respectively, so that there are 
at least 900 spins in the simulation box for each system. 
Each simulation is started with a randomly generated 
spin configuration on the lattice, and 1000 Monte Carlo 
steps were performed to allow the system to equilibrate 
before any measurements were taken. For a system of N 
spins, one Monte Carlo step is equivalent to N Metropo- 
lis loops. (So that in each step, each spin in the system 
is chosen once on average by the algorithm.) After the 
initial equilibration steps, one measurement of energy is 
recorded for each of the 9000 subsequent steps. 10 in- 
dependent simulations are carried out for each extended 
lattice and the results are averaged. 

For the small networks, we carry out 1000 equilibration 
steps followed by 9000 measurement steps as well. How- 
ever, we average the result from 100 sets of Monte Carlo 
simulations so that our results are statistically signifi- 
cant. Because each of the 100 Monte Carlo simulations 
are independent of one another, it is possible to carry out 
the 100 sets of Monte Carlo simulations simultaneously 
with a single processor, which significantly reduces the 
computer time required for the simulation. 



In this paper, we are interested in calculating the en- 
ergy and the specific heat of our systems. The energy can 
be calculated directly by finding the ensemble average of 
the Hamiltonian, which is just the averaged energy for 
all the microstates sampled: 



E = (n^) 



(5) 



where the subscript /i denotes a sampled microstate. 

On the other hand, the heat capacity of the system can 
be calculated directly from the energy by applying the 
fiuctuation dissipation theorem on the system. ^^ As the 
heat capacity is directly proportional to the total num- 
ber of bonds in the system, we have to divide the heat 
capacity by the total number of bonds to allow a fair 
comparison between the various models. In our calcula- 
tion, we define the specific heat, c, as the heat capacity 
per nearest-neighbor bond: 

where N is the number of spins in the system and z is 
the coordination number for each spin. The factor of 2 
corrects for double counting. 

Frustrated Ising systems that have a uniform interac- 
tion constant across the system undergo an "excitation 
transition" as the temperature increases from zero. Be- 
low the temperature of the excitation transition, the sys- 
tem is in one of its ground states. At the excitation 
transition, the excited states become accessible to the 
system. The temperature of the excitation transition is 
largely dependent on the geometry of the Ising system, 
and we have found that the excitation transition of an 
extended frustrated Ising system can be modeled using 
our small networks. The behavior of the excitation tran- 
sition can be best studied using the specific heat as a 
function of temperature, in which the excitation transi- 
tion manifests as a round and broad peak. In addition, 
two systems with the same specific heat vs. temperature 
profile have the same energy profile as well. In view of 
this, we focus on using small networks to approximate 
the specific heat vs. kT profiles of the extended systems. 



B. Analytical Calculation 

As our small networks have a very limited number of 
spins, their specific heat may also be calculated using 
the canonical ensemble and considering all possible mi- 
crostates for the system. In this formulation, the parti- 
tion function is given by: 



-^E^ 



all states 



(7) 



where the subscript ja designates each microstate and 
/3 = 1/kT. The specific heat can be calculated from 
the partition function using the following expression!^ 



Nz 



kf3' 



, S^ In Z 

df3^ 



For very small systems of less than 25 spins, it is easy to 
numerically evaluate the partition function using Eq. ^ 
and then calculate the specific heat using Eq. ([s]). This 
provides us with an alternative way to calculate the spe- 
cific heat of our small networks, and this exact method 
requires much shorter computation time than a Monte 
Carlo simulation. 



C. The Deviation Index 

In this work, we consider a number of small networks 
as approximations to the extended lattices under study. 
Some are better approximations to the corresponding ex- 
tended lattice, while others are not as good. To evaluate 
how good a small-network approximation is, we define a 
deviation index, D, which is a generalization of the coef- 
ficient of determination, i?^, in statistics. ^^ Usually, R^ 
is used to determine the accuracy of a curve fit and has 
values between and 1. The closer the value of R^ to 1, 
the better the approximation. However, in our case, the 
value of D may go below zero when the two curves are 
very far apart. 

Our simulations are performed for the same set of 
kT/J values for each system. For the nth value of kT/J, 
we obtain the specific heat values Cext,n and Csn,n for the 
extended system and the small network, respectively. We 
then calculate the deviation index D using the following 
equations: 



^^tot — / ^ (<^ext,n Cextj 

n 



and 
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SS.. 



ss, 



tot 



(9) 
(10) 

(11) 



(8) 



where / is the mean value of fn. 

To allow comparisons over a range of values of kT / J ^ 
we carry out our simulation for 150 values of kT/J 
equally spaced in logarithmic scale between kT/J = 0.1 
and kT/J = 100. 



IV. COMPARISON BETWEEN SMALL 
NETWORKS AND EXTENDED SYSTEMS 

Geometrically frustrated systems are notable for their 
short correlation lengths. As suggested by Table [T| 
the spin-spin correlation length ranges from extremely 
small for the triangular Kagome lattice to power-law sup- 
pressed for the triangular lattice with the Kagome lattice 
intermediate with exponentially-suppressed correlations. 
It is this absence of long-range correlations that make 
it possible to approximate the collective behavior of an 
infinite system by a small number of spins connected in 



an appropriate network. In a suitable small network, the 
local interactions are preserved but the long-range corre- 
lations are neglected. Accordingly, we expect our small 
network approximation method to succeed for a partic- 
ular extended system precisely to the degree that that 
system has short-ranged correlations. 



Lattice 



Spin-Spin Correlation 



Triangular Lattice 

Kagome Lattice 

Triangular Kagome Lattice 



C(r/ro)"'/^ with \C\ 
e-"/^ with ^ = 3.3r, 
for r/rbb> l f\ 



^ Reference^^ 
b Reference29i3ni 

^ ReferencJ^ 

TABLE L Spin-spin correlations of three antiferromagnetic 
lattices at zero temperature and zero field, where ro is the 
lattice constant for the triangular and Kagome lattcies and 
Thh is the distance between h sites in the triangular Kagome 
lattice. 



A. The Triangular Kagome Lattice 

We begin by developing a small-network approxima- 
tion for the frustrated triangular Kagome lattice (TKL) 
which has the shortest-range correlations of the three sys- 
tems we consider. Unlike the triangular and the Kagome 
lattices, the TKL does not require all bonds to have J > 
to be frustrated. As long as Jaa > \Jah\-, the system is 
frustrated, remaining disordered down to zero tempera- 
ture,^ regardless of the sign of Jab- In addition, earlier 
work by Loh et. al. has shown that the sign of J ah does not 
affect the energy and the specific heat of the system, due 
to a one-to-one correspondence between the microstates 
for two systems with the same Jaa but with Jah differing 
by a sign.^^ Due to the two kinds of bonds in the TKL, 
there are in general two excitation transitions for the sys- 
tem, appearing as two broad peaks in the specific heat 
vs. kT plot. By seeking to reproduce the interactions and 
the connectivities of the lattice sites in the TKL, we ar- 
rived at the "triangular drying-rack network" , which we 
now show is a very good approximation to the frustrated 
TKL. 

A schematic diagram of our small-network approxima- 
tion to the TKL is shown in Fig. ^a). This structure 
has 9 spins and 18 bonds in total. Two a-trimers are lo- 
cated on top and at the bottom of the "drying-rack", 
and three 6-spins are located on the three sides, pro- 
viding the connection between the top and the bottom 
trimers. Fig. [2|b) displays the connectivity between the 
spins in the drying-rack, showing six aa-bonds and twelve 
a6-bonds in the structure. The ratio between the number 
of aa-bonds and a6-bonds in this structure is the same 
as that in the extended two-dimensional TKL. 

In each set of plots in Fig. [3) we compare the specific 
heat vs. kT profiles for the triangular drying-rack net- 
work and the Ising model on the extended TKL with 





(b) 



FIG. 2. (a) Structure of the triangular drying-rack network; 
(b) Connectivity between spins in the triangular drying-rack 
network 



the same Jaa and Jah values. Although the triangular 
drying-rack network has only 9 spins, its specific heat 
profile matches that of the extended TKL almost exactly. 
Since the specific heat profile of the full TKL can be 
modeled by such a limited number of spins in our trian- 
gular drying-rack network, it suggests that it is possible 
to study frustration by focusing on the local interactions 
among a small number of spins that are connected in an 
appropriate network. 



B. The Kagome Lattice 

In this section, we compare the specific heat profiles of 
particular small Ising networks with that of the frustrated 
Ising system on the Kagome lattice. In the Kagome Ising 
system, each spin has four nearest neighbors that are 
arranged in a bowtie network around the spin, and all 
spins on the lattice are geometrically equivalent to one 
another. In view of these geometrical features of the 
Kagome lattice, we consider Ising networks represented 
by two quasiregular polyhedra: the cuboctahedron and 
the icosidodecahedron. The structures of the cuboctahe- 
dron and the icosidodecahedron are shown in Fig. llFa) 
and (b) respectively. We treat each vertex on the polyhe- 
dra as a spin position and each edge as a bond between 
the neighboring spins. Here, the cuboctahedron has 12 
spins and 24 bonds, while the icosidodecahedron has 30 
spins and 30 bonds. Both small networks are similar 
to the Kagome lattice in the sense that each spin and 
its four neighbors form a bowtie network and all spins 
are geometrically equivalent to one another within the 
structure. However, while the Kagome lattice consists of 
triangular and hexagonal plaquettes, the cuboctahedron 
and the icosadodecahedron differ by containing square or 
pentagonal instead of hexagonal faces. 

The specific heat profiles for the Ising systems on 
the Kagome lattice, the cuboctahedron network and the 
icosidodecahderon network are shown in Fig. [5] Compar- 
isons show that the icosidodecahedron's network, which 
has only 30 spins, can serve as a very good approxima- 
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FIG. 3. Specific heat vs. kT plots for tfie triangular drying 
rack network and the extended two-dimensional TKL with 

(a) Jaa/ I Jab| = 5 and (b) Jaa/ \Jab\ = 9. 





tion to the specific heat profile of the frustrated Kagome 
lattice at all kT. On the other hand, the cubotahedron 
network is a less accurate approximation, which produces 
a slightly lower peak than the extended Kagome lattice 
in the specific heat vs. kT plot. However, the cubota- 
hedron network has the advantage that it has only 12 
spins, which makes it very easy to analytically evaluate 
the specific heat exactly using the partition function. In 
Fig. [5) the result for the cuboctahedron network is ob- 
tained using the analytical method, while the result for 
the icosidodecahedron network is obtained from Monte 
Carlo simulation. 



0.06 



0.05 



0.04 



0.03 



0.02 



0.01 



^ extended Kagome 

cuboctahedron 

icosidodecahedron 




10^ 



kT/J 



FIG. 5. The specific heat vs. kT plot for frustrated Ising 
systems on the Kagome lattice, the cuboctahedron network 
and the icosidodecahedron network 

Both the cuboctahedron and the icosidodecahedron 
belong to the class of quasiregular polyhedra, as all the 
vertices are geometrically equivalent to one another on 
each polyhedron. The geometrical equivalence of all spins 
is similarly true for the Kagome lattice. However, we 
show next that it is possible to construct a small net- 
work that well- approximates the Kagome lattice without 
requiring geometrical equivalence between the sites. Con- 
sider two other networks, which we name the "5-bowtie 
network" and the "6-bowtie network", which are com- 
posed of five or six bowties wrapped around in a cir- 
cle. The structures of these two networks are shown in 
Fig. [6J In these two networks, not all spins are geometri- 
cally equivalent to one another. The specific heat vs. kT 
curves for these bowtie networks in Fig. [7| show that they 
are also good approximations to the Kagome lattice. 



C. The Triangular Lattice 



FIG. 4. Structure of Ising systems on (a) a cuboctohedron, 
and (b) an icosidodecahedron. 



The Ising model on the triangular lattice is composed 
only of triangular plaquettes, and each spin is connected 
to six neighboring spins. The best small- network approx- 
imation to the triangular lattice would be obtained from 





FIG. 6. Structure of Ising systems on (a) the 5-bowtie net- 
work, and (b) the 6-bowtie network. 
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mat ion to the extended triangular lattice. However, as we 
consider small networks from tetrahedron to octahedron 
to icosahedron, the approximation gets better. Even 
though there is no network (that we have considered) 
that gives an excellent fit, the approximations given by 
the octahedron and the icosahedron networks are pretty 
good given that they only have 6 and 12 spins respec- 
tively and can therefore be exactly solved with minimal 
computational power. 
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FIG. 9. Specific heat of Ising networks on three regular poly- 
hedra and the triangular lattice 



FIG. 7. The specific heat vs. kT plot for frustrated Ising 
systems on the Kagome lattice, the 5-bowtie network and the 
6-bowtie network 



a polyhedron on which each face is a triangle and each 
vertex is connected to six neighboring vertices. How- 
ever, it can be shown^^ with standard topological meth- 
ods that no such polyhedron exists. Instead, we consider 
networks formed from existing regular triangulated poly- 
hedra, namely the tetrahedron, the octahedron and the 
icosahedron. The geometric structures of these polyhe- 
dra are shown in Fig. [8] 





(a) Tetrahedron (b) Octahedron (c)Icosahedron 

FIG. 8. Ising systems on regular polyhedra 

The specific heat profiles of Ising networks on these 
polyhedra are shown in Fig. |9] The results show that 
none of these networks gives a particularly great approxi- 



The regular polyhedra presented above do not have 
the same number of neighbors for each spin as the tri- 
angular lattice. To find out whether a better approxi- 
mation can be achieved by a structure in which every 
spin has six neighbors, we consider an "icosahedron+2" 
network, constructed by adding two spins to the icosahe- 



dron network as shown in Fig. 10 The resulting network 
has 14 spins in total and every spin in the structure has 
six neighbors. In Fig. [TT] the specific heat profile for 
the "icosahedron+2" network shows that its behavior is 
very far from that of the triangular lattice, although both 
structures are made up of triangular plaquettes and have 
6 neighbors to each spin. This deviation is due to the ar- 
rangement of the triangular plaquettes: while every bond 
is shared by two triangular plaquettes in the triangular 
lattice, some bonds in the "icosahedron+2" network are 
shared by three or five triangular plaquettes. This re- 
sults in a higher density of triangular plaquettes and a 
corresponding increase in frustration. 

To check whether a small network has to be made of 
triangular plaquettes in order to be a good approxima- 
tion to the triangular lattice, we have also considered a 
network on the snub cube, whose structure is shown in 
Fig. 



12 



This structure has 24 vertices and 60 bonds. 
Each vertex in the structure has five neighboring ver- 
tices (as in the icosahedron), but its surface consists of 
squares and triangles. The specific heat of the frustrated 



Ising network on the snub cube is shown in Fig. 13 to- 




FIG. 10. Ising system on the "icosahedron+2" network, an 
icosahedron with two added spins such that every spin has six 
neighboring spins 



0.25 



o extended triangular 
icosahedron+2 model 



0.15 




0.05 



FIG. 11. Specific heat of Ising networks on an icosahedron 
with two added spins and the triangular lattice 



gether with that of the triangular lattice and the icosa- 
hedron network. We find that the snub cube network is a 
much worse approximation to the triangular lattice than 
the icosahedron network, which also has five neighboring 
spins to each spin but is made up of triangular faces only. 
This result suggests that a good approximation-network 
is obtained for the triangular lattice only if the small 
structure is made up exclusively of triangular plaquettes. 



D. Summary and Comparison of Small-Network 
Approximations 

In Table [ll| we tabulate the deviation index D for the 
various small Ising networks considered in the previous 
sections, calculated according to Eq. (11). The value of 




FIG. 12. An Ising system on the snub cube 



extended triangular/6 
snub cube/5 
icosahedron/5 




D provides us with an indicator of the how well each 
small network approximates its corresponding extended 



kT/J 



FIG. 13. Specific heat of Ising networks on the snub cube, 
icosahedron and the triangular lattice 



lattice. A small network can be considered a very good 
approximation \i D > 0.9. For the TKL, we calculate 
D for systems with Jaa/\Jah\ = 9. The values of D sug- 
gests that the triangular Kagome lattice is very closely 
approximated by the triangular drying-rack network, and 
that the Kagome lattice is very well approximated by the 
cuboctahedron, the icosidodecahedron and the bowtie- 
networks. Furthermore, the icosahedron network, which 
has only 12 spins, serves as a reasonably good approxima- 
tion to the triangular lattice. Drawing on our results for 
the specific cases of the triangular, Kagome, and trian- 
gular Kagome lattices, in the next section we summarize 
general criteria for constructing small networks that are 
good approximations to extended frustrated systems. 



V. GENERAL CRITERIA FOR 

CONSTRUCTING SMALL-NETWORK 

APPROXIMATIONS 

In the previous section, we presented comparisons be- 
tween a variety of small networks for each of our proto- 
type two-dimensional frustrated Ising systems. A natural 
question would be: are there any general rules for con- 
structing a small network so that it well- approximates 
the energetics of a particular extended two-dimensional 



extended lattice 


small system 


N 


D 


TKL 


triangular drying rack 


9 


1.0000 


Kagome lat. 


cuboctahedron 


12 


0.9909 




icosidodecahedron 


30 


0.9988 




5-bowtie 


15 


0.9986 




6-bowtie 


18 


0.9966 


triangular lat. 


tetrahedron 


4 


0.4428 




octahedron 


6 


0.8642 




icosahedron 


12 


0.9267 




icosahedron+2 


14 


-1.5526 




snub cube 


24 


-0.4429 




PBC on 4 X 4 unit cells 


16 


0.3416 



TABLE II. Values of coefficient of Determination D for vari- 
ous small networks 



lattice? Based on our results, we suggest the following 
general criteria for constructing small networks to model 
extended two-dimensional lattices: 

1. The local lattice structure around each spin 
in the small network should resemble that in the 
extended lattice. 

By local lattice structure, we mean the arrangement 
of triangular plaquettes around each spin in the net- 
work. For example, each spin in the Kagome lattice is a 
shared vertex of two disconnected triangular plaquettes, 
as shown in Fig. 14]^a). We observe that the better ap- 
proximations to the frustrated kagome lattice, such as the 
icosidodecahedron, the cuboctahedron and the n-bowtie 
networks, have this same local structure around each of 
their spins. Similarly, the triangular drying-rack network 
gives an excellent approximation to the frustrated TKL, 
and the spins on the two structures have identical local 
structures. In both the "drying-rack" network and the 
TKL, the 6- spins have the same local structure shown in 
Fig. 14 ^b) and the a-spins have the same local structure 
shown in Fig. 14 'c). 



V V 

<> <) a(> Oa 



(a) 



(b) 




FIG. 14. Local lattice structures around (a) each spin in 
the cuboctahedron, the icosidodecahedron, the A/'-bowtie net- 
works and the Kagome lattice, (b) each 6-spin in the triangu- 
lar drying-rack network and the TKL, and (c) each a-spin in 
the triangular drying-rack network and the TKL. The orange 
spheres denote the center spins and the grey spheres denote 
the neighboring spins. 

For the triangular lattice, we are not able to find a 
small network that has exactly the same local structure 



as the extended lattice, due to the constraints of topol- 
ogy. However, we observe that a better approximation 
is obtained as the local structure of the small network 
approaches that of the triangular lattice. The triangu- 
lar lattice has a local structure of six connected triangu- 
lar plaquettes around each spin as shown in Fig. [l5|a). 
Among the polyhedra, the local structure of the icosahe- 
dron, with five connected triangular plaquettes around 
each spin as shown in Fig. 15 ^b), resembles the triangu- 



lar lattice to the greatest extent. As a result, the network 
on the icosahedron gives the best approximation to the 
frustrated triangular lattice. On the other hand, the snub 
cube network has a local structure shown in Fig. [l5|c) , 
in which not all the triangular plaquettes are connected. 
As the local structure for the snub cube differs sharply 
from that of the triangular lattice, the behavior of the 
snub cube network is very different from the triangular 
lattice. 






(a) 



(b) 



(c) 



FIG. 15. Local structures around each spin in (a) the trian- 
gular lattice, (b) the icosahedron and (c) the snub cube. The 
orange spheres denotes the center spins and the grey spheres 
denote the neighboring spins. 

We note that the importance of the local arrangement 
of triangular plaquettes results from the fact that the tri- 
angular plaquettes are the main contributors to the frus- 
trated behavior in Ising systems on the triangular lattice, 
the Kagome lattice, and the triangular Kagome lattice. 
The triangular plaquettes give rise to a higher density 
of frustrated bonds than other (2n + l)-gons. If a geo- 
metrically frustrated lattice is made up of pentagons and 
other higher-sided polygons, then the local arrangement 
of the pentagons would be the most important factor to 
consider. 

2. Each bond in the small network should be 
a shared edge of the same number of triangular 
plaquettes as that in the extended lattice. 

In the lattices that we have considered, every edge is 
shared by no more than two triangles due to the con- 
straints of two-dimensionality. When constructing small 
networks for these two-dimensional lattices, we find it 
is extremely important to make sure that every bond is 
shared by the same number of triangles as in the corre- 
sponding extended lattice. For example, in the Kagome 
lattice, every bond belongs to only one triangular plaque- 
tte. In the small networks which give good approxima- 
tions for the Kagome lattice (including the icosidodec- 
ahedron, the cuboctahedron and the bowtie networks), 
every bond also belongs to only one triangular plaque- 
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tte. The same holds true for the triangular drying-rack 
network model for the TKL. 

On the other hand, if this condition is not satisfied, the 
behavior of the small network deviates sharply from that 
of the the extended system. For example, even though 
the number of nearest neighbors are the same in both the 
triangular lattice and the "icosahedron+2" network, the 
later gives a D value as poor as —1.5. The chief reason is 
that there are bonds that are shared by three or five tri- 
angles, which leads to extra frustration due to the excess 
connectivity among those triangular plaquettes. On the 
other hand, in the snub cube network, there are bonds 
that belong to one triangle only, and this also causes the 
behavior of this small network to deviate sharply from 
that of the triangular lattice. 

While the first criterion on the similarity of local struc- 
tures should be satisfied as closely as possible, the second 
criterion that the bonds must be shared by the same num- 
ber of triangles must be satisfied for all small networks. 
Based on our observations, the behavior of a frustrated 
lattice is heavily dependent on the types of bonds in the 
network. Whether a bond is shared by one, two, or more 
triangular plaquettes makes a drastic difference, as the 
connectivity of triangular plaquettes through the bonds 
has a significant effect on the density of states and the 
number of accessible states in the system. For example, 
in the ground states of the snub cube, the square pla- 
quettes are not frustrated, and thus the snub cube is less 
frustrated than the icosahedron, even though both struc- 
tures have the same number of neighbors for each spin. 
The snub cube thereby has a lower energy per spin in the 
ground state, resulting in a much higher peak in the spe- 
cific heat than the icosahedron (whose behavior is close 
to that of the frustrated triangular lattice) . 



VI. SMALL NETWORKS VS. PERIODIC 
BOUNDARY CONDITIONS 



For extended two-dimensional systems whose spin-spin 
correlation decays with distance r, it is difficult to know 
a priori the minumum size required for a patch of the lat- 
tice with PBC to well- approximate an extended system. 
We now show that (for a comparable number of spins), a 
network constructed by applying PBC to a small section 
of the extended lattice gives a much worse approxima- 
tion than the small networks we have developed. For the 
triangular lattice, applying PBC on 4 x 4 unit cells (16 
spins) is insufficient to produce a good approximation to 
the extended triangular lattice. Fig. p!6] shows the specific 
heat profile obtained from simulations on 4 x 4 unit cells 
of the triangular lattice with PBC, as well as the specific 
heat profiles for the octahedron network, the icosahedron 
network, and the triangular lattice for comparison. Even 
though 4x4 unit cells has more spins than the octahe- 
dron and the icosahedron, its fit of the triangular lattice 
has a D value of 0.3416 (Table [TT|) , much worse than the 
icosahedron and the octahedron. The reason behind this 
huge deviation is that by applying PBC, one introduces 
additional triangular plaquettes. 



extended triangular/6 
octahedron/4 
icosahedron/5 
4x4 unit cells/6 




kT/J 



In Section [V| we presented general guidelines for con- 
structing small networks to approximate the thermody- 
namics of extended two-dimensional lattices. Our de- 
gree of success is directly related to the short correlation 
lengths in these systems, which allow local structure to 
dominate the thermodynamics. Given this, the reader 
may be wondering whether just as good of an approxi- 
mation can be obtained by taking a small piece of the 
two-dimensional lattice and applying periodic boundary 
conditions (PBC). If the correlation length is zero outside 
of a unit cell (as it is for the triangular Kagome lattice), 
applying periodic boundary conditions on one unit cell 
would be expected to give a very good approximation. 
In fact, the network obtained from a unit cell of the tri- 
angular Kagome lattice with PBC is identical to what 
we have called the triangular drying-rack network. How- 
ever, applying PBC to the Kagome and triangular lat- 
tices, which have longer correlation lengths, is met with 
much lower success, as we show next. 



FIG. 16. The specific heat profile obtained by applying peri- 
odic boundary condition on 4 x 4 unit cells of the triangular 
lattice compared with that for the extended triangular lattice, 
the octahedron and the icosahedron 

We conclude that, compared to applying PBC to a 
finite piece of the lattice, the method of small-network 
approximation that we have developed in Section |V| is 
generally more accurate, for the same number of spins. 
In some cases, one of our small networks has the same 
structure as a few unit cells of the extended lattice un- 
der PBC. In addition to the case of the TKL and the 
triangular-drying-rack-network, for the Kagome lattice, 
the network of the cuboctahedron turns out to be the 
same as applying PBC on 2 x 2 unit cells of the Kagome 
lattice. However, the variety of structures covered by our 
small-network approximation is much larger than what 
can be obtained by applying PBC and in the case of the 
triangular lattice, the small-network approach allows us 
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to do much better than applying PBC. 



VII. CONCLUSION 

In this work, we have examined the correspondence 
between geometricahy frustrated Ising systems on par- 
ticular small spin networks and three ordinary extended 
two-dimensional systems. The small correlation lengths 
in these frustrated Ising systems makes it possible for 
a remarkably small network of spins to give a good ap- 
proximation for the specific heat of the corresponding 
extended two-dimensional lattice. However, the correla- 
tion length poses a fundamental lower limit on the size 
of the small network required to obtain a good approx- 
imation, and we find decreasing success for this small- 
network approach as the correlation length of the ex- 
tended lattice system increases. We suggest that using 
suitably-designed small networks is a good way to ob- 
tain a first approximation for the properties of large and 
complicated geometrically frustrated systems. 

The behavior of frustrated magnetic systems are of- 
ten difficult to study, particularly beyond the simplest 
approximations. For many physical systems, Heisenberg 
spins are required to accurately model the quantum prop- 



erties of the systems. Even within the classical Ising 
approximation, computationally-intensive Monte Carlo 
simulations on large systems are often used because there 
are no general analytic approximations available for frus- 
trated systems. In this work, we have presented an ap- 
proach that may provide a general way to approximate 
frustrated systems with extremely small computational 
expense. The small networks discussed here could be em- 
bellished to resemble more realistic models, for instance 
by using Heisenberg spins instead of Ising ones, while re- 
maining computationally feasible. In addition, the idea 
of small- network approximations could lead to new renor- 
malization or numerical procedures. In this work, we 
have focused on ordinary two-dimensional systems that 
are frustrated due to their triangular plaquettes. How- 
ever, similar methods could be applied to other kinds of 
geometrically frustrated lattices with various geometri- 
cal structures in both two- and three-dimensions, such 
as spin ice. 
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